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The gamma-ray binary system PSR B1259-63 has recently passed through periastron 
and has been of particular interest as it was observed by Fermi near the December 2010 
periastron passage. The system has been detected at very high energies with H.E.S.S. 
The most probable production mechanism is inverse Compton scattering between target 
photons from the optical companion and disc, and relativistic electrons in the pulsar 
wind. We present results of a full anisotropic inverse Compton scattering model of the 
system, taking into account the IR excess from the extended circumstellar disc around 
the optical companion. 
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1. Introduction 

An important Be X-ray Pulsar binary (Be-XPB) example is the TeV gamma-ray 
system PSR B1259-63, which consists of a 48 ms pulsar orbiting the Be star LS 
2883P The orbit is eccentric (e « 0.87) and unpulsed TeV gamma-ray emission 
has been detected by the High Energy Stereoscopic System (H.E.S.S.) close to pre- 
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vious periastron passages. 1 ^ 1 ^ The circumstellar disc is mis-aligned to the orbital 
plane and the pulsar passes through it twice, before and after periastron. Previ- 
ous observations have also detected unpulsed radio emission near periastronpJ and 
X-ray emission across the whole orbit!^ The system is powered by the spin-down 
luminosity of the pulsar and the unpulsed radiation is believed to originate from a 
stand-off shock-front that forms between the pulsar and stellar winds. The X-ray 
emission is probably synchrotron radiation and gamma-rays are produced through 
inverse Compton (IC) scattering between target photons from the Be star-disc sys- 
tem and relativistic electrons (7 ~ 10 6 )Pln this study the anisotropic IC scattering 
in PSR B1259-63/LS 2883 has been investigated, focusing on the contribution of 
the infrared (IR) excess from the circumstellar disc to the IC scattering rate. The 
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additional effec ts o f photon absorption through pair-productiorPI or changes in the 
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shock condition' ° have not been considered here. This study only presents the ef- 
fect of the IR excess to investigate whether this a viable source of additional photons 
and whether this will have a significant influence on the gamma-ray production. 

2. Modelling 

In order to model the anisotropic IC spectrum from the circumstellar disc in 
PSR B1259-63 it is necessary to calculate the total number of scatterings, tak- 
ing into account the anisotropic scattering rate and the photon contribution from 
the star and the circumstellar disc. The modelling presented here is dependent on 
the IR excess and curve of growth modelling calculated in Ref. [5] which assumed the 
stellar parameters presented by Ref. [TUJ For this reason a stellar mass and radius of 
IOMq and QRq, and a pulsar mass of IAMq have been assumed and not the newer 
system parameters presented by Ref. [TTJ 

2.1. Anisotropic Scattering 

The anisotropic scattering rat JEH31 neec j s t take into account the scattering angle, 
#0) between the direction of the observer and the incoming photons over the full 
solid angle of the star and disc. The total anisotropic scattering rate is given bjD31 



where n p h is the number of photons, n e is the electron distribution, dN/dtdeo is the 
anisotropic scattering rate dependent on the scattering angle 9q, and cos 9 is the 
angle to the normal of the photon source (star or disc). The star is assumed to have 
a black-body photon distribution (T» = 25 000 K), while the photon distribution 
from the disc is calculated using the curve of growth method. 

The electron spectrum is assumed to have the form n c (-y) — K c ~/~ p , where K e = 
1 is used in this study as only the resulting change due to the IR excess is considered. 
The synchrotron radiative cooling time is t sync ~ 770 (7/10 6 ) _1 (.B/l G) -2 seconds, 
which is greater than the light crossing time of the disc. For this reason only a 
constant electron spectrum has been considered, appropriate for adiabatically cooled 
electrons!^ Therefore, this study considers only whether the IR excess from the 
large circumstellar disc significantly influences the IC scattering. Additional effects 
due to radiative cooling have not been considered here. 

2.2. Disc photon spectrum - curve of growth method 

The intensity from the circumstellar disc depends on the free-free optical depth 
through the disc, from which the photon number d ensi ty is calculated. The optical 
depth is determined by the curve of growth method^ and the pulsar-disc orienta- 
tion. It is assumed that the circumstellar disc has a half-opening angle #disc, extends 
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Fig. 1. (Left) Scattering rate for the disc and star contributions in PSR B 1259-63 at periastron 
for a disc orientated 90° and 45° to the orbital plane, scaled to the peak stellar contribution. 
(Right) The corresponding light curve of the high energy contribution (~ 0.1 — 100 GeV) of the 
star and disc around periastron considering only adiabatically cooled electrons. 



to a radius i?disc and has a power-law density profile that decreases with distance 
as p(r) = po (r/R*)~ n . 



3. Results 

Two orientations of the disc have been considered; the first with a disc tilted 90° to 
the orbital plane with the normal lying along the semi-major axis (Tdi sc = 12 500 K 
& #disc = 5°), the second titled 45° to the orbital plane with th e norm al rotated 
19° from the semi-major axis (T disc = 15 000 K & 6> disc = 0.7°).EHS These disc 
temperatures and small half-openin g angle s are consistent with current modelling 
of the circumstellar discs of Be starsllSllSJ The disc is assumed to have a radius of 
50i?*, consistent with the binary separation at the time of the binary eclipsePSEl 
The parameters for the curve of growth method were taken from the fit to previous 
optical and IR observations^ and the electron spectrum is modelled using an index 
of 2.4 between 7 = 2.54 x 10 5 and 2.54 x 10 7 (adopted from Ref . fl"5|) . Fig.Q]shows the 
relative scattering rate at periastron scaled to the peak of the stellar contribution 
for both disc orientations and the corresponding light curve at high energies (~ 
0.1 — 100 GeV). In both cases the contribution of the disc occurs towards high 
energies (as opposed to very high energies) where it results in an increase in the 
scattering rate by a factor ~ 2. The light curves show a peak in the disc contribution 
near periastron and less contribution near the disc crossing where the IC scattering 
is dominated by the scattering of stellar photons. The rotation of the disc from the 
semi-major axis results in a slightly higher flux before periastron and the shifting 
of the peak away from periastron, while the post-periastron flux is lower because 
of unfavourable scattering angles. In both cases the flux increases after the disc 
crossing as head-head scattering readily occurs in this region. 
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4. Conclusion 

The IC scattering from the circumstellar disc in PSR B1259-63 is shown to increase 
the flux at high energies by a factor of ~ 2 during periastron. Due to the large 
size of the disc and the increase in the photon density at IR wavelength (increasing 
by ~ 1000 at mid-IR wavelengths it was possible that the increased solid angle 
could have a major influence. However, the flux from the circumstellar disc drops 
off rapidly with radial distance and the low flux from the disc at the point of the 
disc crossing mitigates the influence of the disc. This approximation, considering 
only adiabatically cooled electrons, shows that the influence of the IR excess does 
not dramatically increase the scattering rate during the disc crossing. In particular, 
this simplified model does not explain the flare observed by Fermi after the second 
disc crossing It should be noted that a simplified approximation was used to 
calculate the photon spectrum from the disc and a more detailed model may result 
in a higher photon density at the outer edges of the disc. Additional shock heating 
effects may also increase the brightness of the circumstellar disc during the pulsar 
passage. Such an effect could increase the scattering rate near the disc cr ossi ng and 
IC scattering of pre-shock 7 = 10 4 electrons may explain the Fermi flareP^ 
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